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1. Introduction 

Often in Physics at the nanoscale one deals with systems which can be regarded as 
having support in a very small neighborhood C of a one dimensional graph V. 
The parameter £ is a measure of the distance of dV^ from F as compared to the length 
of a typical edge of the graph. 

The dynamics of such systems is described by a propagation equation, which can 
be the wave equation in the case of optical wires or the Schrodingcr equation in the case 
of conduction electrons in a macromolecule or in the case of conducting nano-devices. 
In the case of optical wires it is reasonable to choose Neumann (reflecting) boundary 
conditions at the boundary dV^ . In the case of macromolecules or nano-devices on the 
other hand one can consider that the system is confined in the region F^ by very strong 
forces at the atomic scale, which provide a potential having the form of very deep and 
narrow valleys. In this case the use of Dirichlet boundary conditions can be regarded 
as a reasonable approximation. 

The inital value problem is well posed for both boundary conditions, but in general 
one is not able to give the solution in an explicit form. It is then useful to search 
for a dynamical system which is explicitly solvable and at the same time provides an 
approximation to the physical one when e is very small. In this case it is important 
to give an estimate of the error one makes in this approximation for the observables of 
interest. 

One can expect that the limit dynamical system be described by a wave equation 
or respectively a Schrodinger equation on the graph. But if this is the case, in order 
to have a well posed problem one must specify boundary conditions at the vertices. So 
the problem becomes: which are the boundary conditions (if any) that one must choose 
at the vertices of the (metric) graph, in order to have an evolution equation on the 
graph with solutions which give an approximation (to order e", for some a > 0) for 
the expectation values of the relevant physical observables? One may reasonably expect 
that the answer depends on the shape of F*^ in the vicinity of the vertices. 

In the case of Neumann boundary conditions the answer is known (Exner & 
Post 2005) for the case of the Schrodinger equation, at least for initial data with not too 
large energy. Indeed in this case the initial wavefunction can be chosen smooth (e. g. 
in W"^'"^) uniformly in e and the solution can be restricted to the graph uniformly when 
£ — )■ 0. As a consequence, one can define a sequence of maps from F^ to F, and prove 
that the trace on F of the resolvent of the Schrodinger operator on F^ converges to the 
resolvent of a Schrodinger operator on F with Kirchhoff (coupling) vertex conditions 
(continuity of the solution, zero sum of the outward derivatives). The situation is 
entirely different in the case of Dirichlet boundary conditions on dV^: the VT^'^-norm of 
the initial datum increases without bound as £ — )■ 0, and there is no limit trace on F. 
For this reason, the question of the existence of a limit flow had not been answered so 
far (Kuchment 2008), in spite of the obvious physical and mathematical interest. 

In concrete cases of physical devices, models of limit dynamics have been 
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constructed to fit the experimental data. For instance, in the case of a sharply 
bent conducting device experimental evidence shows that in order to have a good 
approximation one must use Dirichlet (decoupling) boundary conditions at the bend. On 
the other hand, the standard treatment of conducting electrons in aromatic molecules 
(such as graphene or benzene), in which the molecule is represented by a graph, shows 
that results in fair accordance with experiments are obtained if the limit model is 
constructed with conditions at the vertices that are of weighted Kirchhoff type. The 
phenomenological values of the weights are different for different molecules. One can 
consider this as an evidence that the "right" boundary conditions depend on the shape 
of 

Indeed the images seen at the electronic microscope are different for different 
molecules (one takes to be the region in which the density of conducting electrons 
is appreciably different from zero) but in all cases they have the shape of an annulus 
around the nucleus, with tunnels in correspondence of the valence bonds. One interprets 
this structure as due to the combined action of the attraction by the nucleus and the 
exclusion principle that forbids the conduction electrons to occupy the region of the core 
electrons. 

Here we give an answer to the mathematical question described above. We show 
that the generator of the limit dynamics on the graph corresponds to a suitable boundary 
condition at the i-th vertex Vi, which depends on the shape of F^ near Vi through the 
spectral properties of a sequence of auxiliary Schrodinger operators defined on a suitable 
neighbourhood of the vertex Vi. 

In this paper we analyze the problem of approximating the free Schrodinger 
evolution on an e-thin Dirichlet star-shaped waveguide, with semi-infinite cylindrical 
ends of thickness £; in a future publication we shall generalize our results to cover generic 
waveguides and smooth potentials. In this paper we study convergence of solutions and 
(weak) convergence of resolvents, and briefly suggest how to add an external potential 
and study convergence of scattering matrices. 

2. Description of the problem 

In this section, we describe the structure of the waveguide in detail, and we introduce 
some notation. The waveguide F^ is obtained by gluing smoothly to the vertex region 
Jl^ (a compact set in R^) a number of semi-infinite cylindrical ends, the branches of 
the waveguide; the boundary dV^ is smoothf. The cylindrical branches have two- 
dimensional section E^, a compact domain in with smooth boundary, of linear 
dimension proportional to e. Also, the linear dimension of the vertex region is 
proportional to e so that the waveguides {F^} associated to different values of e are 
connected by a similarity transformation (a change in the length scale). With this 

I For illustrative purposes, in the Figures we depict two-dimensional waveguides with edges at some 
points of the boundary. Note that the hypothesis of smooth boundary could be weakened, and the 
dimension of the waveguide could be any d>2. 



Effective dynamics on thin Dirichlet waveguides 



4 



assumption, the geometry of is essentially fixed, and the problem is equivalent to the 
computation of the low-energy effective dynamics in a fixed domain. Note that this is 
not true anymore when we have many vertices at a fixed distance, and we shrink the 
cylinders connecting them. 



Structure of the waveguide 

The branches of the waveguide are labeled by the index j, which runs from 1 to n. Each 
branch is isometric to the cylinder E*^ x [0, +oo); then 

/ n ^ 
r = U IJ S^^ X [0, +oo) 

The Dirichlet Laplacian on has discrete spectrum 

{/^l,/^2---} fll<fll<... 

and the eigenvalues satisfy the scaling relation 

The Dirichlet Laplacian on the entire waveguide has absolutely continuous 
spectrum which coincides with the semi-infinite interval [/xf,+oo), and possibly some 
discrete eigenvalues Af . . . A|, not larger than fif, corresponding to bound states§ with 
eigenfunctions 0f . . . 0|. They also satisfy the scaling rule 

\£ — 

'^i ~ 72 ■ 

It is natural to work with the renormalized Hamiltonian 
= -Ac - ^\ 

acting on L^(r^). In this way, the bound states have negative energy and the absolutely 
continuous spectrum coincides with the positive half-line. 

Let us introduce the crucial concept of mesoscopic region F^^^j, or sometimes simply 
Tint- This is obtained by cutting at distance i from the vertex region fl'^: 

r-;; = ^^^u(Us^^x[o, 

The complement of Tint is referred to as Tout- 

n 

^'oL = [j^x % +oo) . 

It is very convenient to work with the parameter L defined as 

L = £/£ 

§ bound states are known to form e.g. near the bent regions of the waveguides; see for example (Exner 
et al. 1989) 
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instead that with We refer to L as to the rescaled mesoscopic lenght. With this 
notation, the mesoscopic region has tickness e and size Le. 

Remark To identify the points in the j'-th branch of the waveguide, we normally 
use Cartesian coordinates {xj,yj), with values in [0, +00) x S^. We can say that these 
coordinates are adapted to the macroscopic scale. Later it will be convenient to make 
use also of rescaled coordinates {Xj,Yj), adapted to the microscopic scale, ranging in 
[0, +00) X S^, such that 

Xj = sXj, yj = sYj . 

Notice that in these coordinates, the mesoscopic region is characterized by the 
inequalities 

Xj <e, Xj<L. 



Le 



Figure 1. A very simple example of waveguide, the L-shaped waveguide. The junction 
is a square of size s, and there are two semi-infinite rectangular branches attached 
to it. The mesoscopic region T^^^ is white, while F^^fj is shadowed. 



The general strategy The waveguide that we consider is characterized by two 
parameters: the microscopic lenght e and the mesoscopic lenght Le — £. The parameter 
i is chosen so that the inequality 

1 

is satisfied: this motivates our notation, too. In terms of the rescaled mesoscopic lenght 
L, these conditions read 

£<1 L>1 L-£<1. 
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Our strategy is the following: given a wavcfunction if) on F^, solution of the 
Schrodinger equation, we consider its restriction to F^^f^ and F^^^. We prove that ip \ F^^ 
essentially factors into ®Xi plus negligible corrections: \1' is the effective wavcfunction 
on the limit graph F. It is clear that this approach makes sense if £ <C 1, so that ^ is 
defined everywhere except for a small neighbourhood of the vertex. 

Our second step is to understand the behaviuor of ^ near the vertex. To do this, 
we study the restriction of to F^^^. It turns out that for a suitable class of initial 
states (roughly speaking, low-energy states) the behaviour of ^ at the vertex does not 
depend on the initial state, but only on the spectral properties of F^^^^. The estimates 
that constrain the behaviour of ^ are expressed in terms of both the parameters i, e: 
in concrete situations one can choose £ to be a suitable function of £, to optimize the 
error. This is particularly simple for free particles, but we prefer to keep the discussion 
at a general level to prepare the ground for successive generalizations. 




Figure 2. A more complicated waveguide. Three branches of thickness e join at a 
junction 0,"^ . In this case, the junction is the inner part of a circular crown, of radii 
TiE, r-iE- The mesoscopic region is obtained by attaching to the junction a portion of 
the leads of lenght Le. In general, one may consider even more irregular junctions, as 
long as they are obtained by rescaling a fixed compact domain by e. 
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Some definitions 

We introduce here the notion of spectrum of the mesoscopic region, and illustrate some 
results from (Grieser 2008a) that characterize it. Given a generic junction , it is clear 
that the explicit spectrum cannot be computed in terms of the parameters e and L, even 
if some asymptotics may be calculated for specific geometries. Our point of view is that 
we can formulate some hypotheses on this spectrum regarding its qualitative behaviour 
near the continuum treshold, that enable us to distinguish between junctions that lead 
to decoupling or coupling conditions for the limit Quantum Graph. 

Definition 1 (auxiliciry Hamiltonian) On the mesoscopic region T^^^ we define the 
operator 

where the Laplace operator A is defined by Neumann boundary conditions at {£} x 
and Dirichlet boundary conditions on the rest ofdF^^^. H^^^ is the auxiliary Hamiltonian 
on the mesoscopic region. The spectrum of H^^^ is discrete and it will be denoted by 

♦ 

The spectrum of the auxiliary Hamiltonian H^^^ is sometimes referred to as the 
spectrum of the mesoscopic region. This spectrum is analyzed in (Grieser 2008a) in 
great detail, although in a slighly different setting: the domain is rescaled by a factor 
1/e and the eigenvalues are multiplied by accordingly. The dilated mesoscopic region 
has linear dimension L and the branches have tickness 1; sometimes one just says that 
we are working in coordinates such that e = 1. In particular, it is shown that any 
eigenvalue (proper and generalized) of H'^ can be approximated by an eigenvalue of H^^^ 
and vice-versa. For the point spectrum, the difference between A| and A^'^ vanishes 
exponentially in L = l/e: in our notation, it is 0{e~'^e~^). The eigenfunctions are 
shown to converge too, in the proper topology, with the same speed. 

An immediate consequence of the considerations above is that the first k eigenvalues 
of Hl^^ will converge to the energy levels A^ . . . A| corresponding to the bound states of 
r^; the other eigenvalues of Hl^^ will "merge into the continuum". 

For each value of e, we shall divide the eigenvalues that merge into the continuous 
spectrum in two classes, according to qualitative properties of the approach to the 
bottom of the continuous spectrum when e — > and L — )■ oo. If there arc no elements 
in the first class, the limit motion on the graph has Dirichlet (decoupling) boundary 
conditions at the vertex of the graph. The eigenfunctions corresponding to the elements 
of the first class (we call them resonant sequences) provide "bridges" connecting the 
different branches of the graph; if there is only one element in this class, the motion 
on the limit graph is generated by a Laplacian with suitable Kirchhoff-like boundary 
conditions, which depend on the asymptotic form of this connecting eigenf unction. If 
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there are more than one elements, the problem becomes more complex, and will not be 
treated here. 

Recall that the first k eigenvalues A^'^ . . . A^'^ of H^^^ converge to the k bound states 
A^ . . . A| of the waveguide as e — ?> 0; henceforth, we will loosely refer to all these states 
as "bound states", or "localized states". To proceed further, it will be convenient to 
introduce the following 

Definition 2 (spectral gap condition) The family of (mesoscopic) regions T^-f^ 
satisfies the spectral gap condition if there exist /io > and 7 > such that the 
eigenvalues of the non-bound states ipf^^, n > k + 1 satisfy, uniformly in s > 

♦ 

We are also interested in waveguides where the gap condition is not satisfied, i.e. 
there are some non-bound states with eigenvalues that do not satisfy the inequality 
stated above. We do not treat the general case, but limit ourselves to the following 
important special case: 

Definition 3 (resonant sequence condition) The family of (mesoscopic) regions 
r^^f satisfies the resonant sequence condition if the eigenvalues of the non-bound states 
tpn^, n> k-\-2 (hence, all of them but the {k-\- l)-th) satisfy condition (1) uniformly in 
£ > 0, while 

-§-(^)''<Aii,<0 (2) 

for some 7' > 1, > 0. 
♦ 

We call the sequence of states {V^^+i} which characterize the definition above a 
resonant sequence of states. The name we have chosen comes from the fact that if the 
waveguide has a zero-energy resonant sequence (a solution of the Dirichlet problem not 
in L^(r^)), then its restriction to Vl^^ provides a resonant sequence. The last statement 
follows from a variational argument, which bounds the eigenvalue from above, and a 
result discussed in (Grieser 2008a) which bounds the eigenvalue from below. Note that 
(Grieser 2008a) proves that in the self-similar case the convergence is exponentially fast. 
Indeed, including the polynomial bound is useful for generalization to the non-uniform 
case. 

In a self-similar family of waveguides, if a resonance exists for one value of £, it exists 
for all values of e. But a resonant sequence may exist even if there is no zero-energy 
resonance: this will be very useful for generalising the analysis to graphs with many 
vertices, where self-similarity is lost. We expect the resonant state condition to hold 
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when the waveguide has some resemblance with the domain accessible to conduction 
electrons of an aromatic molecule||. 

Remcirk The inequalities of Equations 1, 2 arc written in terms of the 
parameters L so that their application in the forthcoming Theorems is immediate. 
But we stress that if we rewrite the Definitions above for the rescaled mesoscopic region 
of tickness 1 and lenght L ("take e = 1"), for which the spectrum is rescaled by a factor 
then the eigenvalues are functions of the parameter L alone, and Equations 1, 2 
become conditions on the asymptotic behaviour of these eigenvalues, for L large. Thus, 
the information about the limit Quantum Graph (wether the junction gives coupling or 
decoupling gluing conditions) is encoded into the asymptotics of the spectrum of the 
rescaled mesoscopic region. The parameter e will play a role when we discuss the class 
of states (the low-energy states) that can be approximated by effective wavefunctions 
on the graph; at that stage, it will be convenient to choose L as a function of e in order 
to minimize the error terms. 

Towards the construction of examples It would be very interesting to build some 
concrete examples of geometries such that the Definitions above are fulfilled. We give 
some suggestions that come from preliminary computations. 

We expect that the L-shaped waveguide of Figure 1 falls into Definition 2. This 
waveguide has one bound state (Exner et al. 1989) (in our notation. A; = 1), hence we 
need to look at the behaviour of the second eigenvalue X^'^ of the auxiliary Hamiltonian. 
For this particular geometry it should be possible to compute the asymptotic of the 
spectrum for large L, by adaptating the arguments used in the paper (Grieser & 
Jcrison 2007): we expect that in this case the condition of Equation 1 is fulfilled. 
Note that if iIj2^ is antisymmetric with respect to the reflection through the axis of the 
waveguide, this follows immediately from a standard bracketing argument. 

We also expect that the waveguide depicted in Figure 2 falls into Definition 3 
instead, at least for a wise choice of the radii ri,r2 of the junction. The idea is that 
there are no bound states, so that if we are able to tune the parameters in such a way 
that a resonance for the infinite waveguide is present, then the resonance induces a 
resonant sequence by the mechanism described previously. Note that the "hole" is a 
suggestion that comes from Physics: it mimics a repulsive potential, which is given for 
example by the core electrons of the carbon atom that sits at the vertex of a graphene 
layer. 

Clearly, it is more difficult to exihibit examples of guides satisfying this hypothesis. 
It is to be expected that a resonance (or its counterpart, the resonant sequence) 
is not present for a generic waveguide, accordingly to the general wisdom that 
decouphng generically occurs for thin Dirichlet waveguides (Molchanov & Vainberg 
2007) (Molchanov & Vainberg 2008) (Grieser 2008a). 

II The role of a resonant structure has been advocated by B. Pavlov, from a different point of view. 
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3. Multiscale decoupling 

In this section, we want to show that outside the mesoscopic region low-energy 
wavefunctions are confined to the first transverse mode, and the coefficient solves the 
free one-dimensional Schrodinger equation; therefore, they can be identified with one- 
dimensional wavefunctions on F. The next sections are devoted to describe the behaviour 
of the solutions inside the mesoscopic region, that is to say, near the vertex of F. 

Let -0^ belong to the image of the spectral projector P{o,e] associated to H^. For 
practical purposes, we may assume that at time zero -0^ is approximated by a product 
state 

^l)%t = 0) ~ $ (g) x! 

supported in F^^f^. This state has finite energy, and is orthogonal to the bound states 
0f ... 01 of F^: both these properties hold at all times, by unitarity of the time evolution. 
At any time, the restriction of to T^^^ admits a convergent expansion in Fourier modes: 

oo 

out ' * 

m=l 

where is the m-th Dirichlet eigenfunction of —A on E^, and is a vector- valued 
function on R+; its components are ^^j, the index j runs from 1 to n and labels the 
branches of the waveguide. It is useful to think about for any m, as a function on 
the limiting graph F. 

We will prove the following 

Theorem 1 (Multiscale decoupling) Consider the wavefunction tp^ described above, 
and restrict it to F^^j. We can prove that 

oo 
m=2 

for some numerical constant C. Moreover, 

oo oo 

E ll*^ll^^.^{r) <C'-e, E ll*-lli-{r) < ^ ■ e 

for some numerical constant C. 
♦ 

In other words, the components of the wavefunction along the higher transverse 
modes are suppressed in the limit £ — > 0: the only component which survives is "frozen" 
in the first transverse mode Xi- 

Proof 

The first k eigenvalues X\'^ . . . A^'^ of the auxiUary Hamiltonian H^^^ converge to the 
bound states . . . A|, in the sense explained above: the difference is C(e~^e~^) 
(we will loosely refer to them as "bound states" too); the eigenvectors ipl'^ ■ ■ • V'fe^ 
converge uniformly, with all derivatives, to the respective bound states 0f . . . 0| of 
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restricted to F^^^^. It follows that if -0 is orthogonal to . . . (/)|. , its restriction to 
-'^irat ^i^^ almost orthogonal to ^^'^ . . . V'^,' , up to an error that is 0{e^'^e^^). 
In the following we will just say that i/; (restricted to T^^^J is almost orthogonal, 
or quasi-orthogonal, to Vi'^ • • • V'fe^) without writing down the exponentially small 
error: this is because we will obtain estimates that are polynomial in e and ^, so 
that the exponentially small error is easily reabsorbed; and this will help a lot in 
keeping the formulas clean. 

If the Laplacian on T^^^ has no eigenvalues below /xf other than the "bound states" , 
a standard spectral argument implies that 

/ (|vvi'-/xf|vi') >o. 

Moreover, by substituting the Fourier expansion of |~ we easily see that 



■ out 



Hence, both these contributions are positive: from the energy bound on ip^ we know 
that their sum Jp^ (jVV'^P — /if|^''p) is positive and smaller than E, and so each 
one of them (in particular, the integral over Tout) is bounded by E. Notice that if 
the integral over Tint had been large and negative, the integral over Tcmt could have 
been arbitrarily large; this is why we must avoid bound states. 

After this simple, but important observation, the first part of the theorem easily 
follows by substituting the Fourier decomposition of \ T^^mt i^ integral. We 
obtain 



-lit V/^2-/xi7re 
To complete the theorem, note that HiJ^-^^lp < E"^. Repeating the previous 

argument, with simple modifications, one proves that the sum of the norm squared 
of A^*^ is bounded: by interpolation, the sum of the squared norms of V^^, for 
m > 2 is 0{e). This implies that the VF^'^-norm of m > 2, is 0{e^^'^): we recall 
that this implies that the L°°-norm is infinitesimal too, with the same bounds. 

Remcirk Consider the expectation value of an observable, which is either given 
by the multiplication by a smooth function supported in F^^^, or a linear differential 
operator in the longitudinal derivatives d/dxj, j = 1 . . .n with coefficients in the smooth 
functions supported in F^f^. Theorem 1 tells us that the expectation values on ip"^ of these 
obscrvablcs converge for small e to the expectation values on of the corresponding 
observables restricted to F. Consider for instance a function (p on F^ (vanishing on T^^^): 
and say that we want to compute 

By Theorem 1, it is easy to see that apart from an error term vanishing as a suitable 
power of e, these expectation values respectively converge to 
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where 0r is the evaluation of on F. Notice that ||xilP = 1 Vs. 

The last remark is crucial: it allows us to use instead of if)^ to compute the 
expectation values of important observablcs, such as the density and (longitudinal) 
current of particles. Henceforth, we will study the fundamental Fourier component 
\l/f, and systematically neglect the higher-order contributions to '0'^, which do not give 
contribution to the observables of interest. 

It is desirable to determine an effective evolution equation for ^\ in the hmit £ — > 0. 
Notice that defined as 

is measurable, as its time and space derivatives; the L^-norms of \E'f and of V\E'f are 
bounded, hence ^\ is in Vr^'^(F). It is easy to see that satisfies the free (one- 
dimensional) Schrodinger equation, in weak sense. Indeed, for any smooth test function 
</? with support in F^^, and of the form ®Xi the equations of motion read 

and after the integration on the transverse variable the equation reduces to 

This equation is the weak form of the free evolution of along the edges of the graph. 
Since the support of $ is at a distance not smaller than I from the vertex, the equation 
is satisfied outside the vertex. By a compactness argument, there exists a subsequence 
{^i"}ra6N having a weak limit ^ e l^^'^(F), which must satisfy the weak Schrodinger 
equation on the edges of F as well (recall that the initial datum is the same for all 
values of £!). ^ is our limit wavefunction on F: once we have proven that the unitary 
evolution of ^ is uniquely determined, it will be a posteriori clear that any subsequence, 
and hence the whole sequence {^'i} converges to the same Our goal is to use ^, 
instead of , to approximate the expectation values of ip^ . The next section is devoted 
to determining the time evolution of ^, which depends crucially upon the behaviour of 
^ at the vertex of F. 

4. The limit motion on the graph 

In the following, we will show how one can approximate the low-energy dynamics of 
the dynamical system (F^,i7^) with a suitable dynamics Hy on the metric graph F. 
The cornerstone of our approach is the study of the spectrum in the mesoscopic region 
F^^^j. The limit dynamics is uniquely identified in this section; more quantitative results 
(convergence of resolvents and of wave operators) are dealt with in the next sections. 

For simplicity, we only treat the case with no external potential. The method can 
be used also in case there is a continuous potential V{x), x e IR^; we shall denote 
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by Vy its restriction to the graph. We suggest that such perturbations can be treated 
by approximating the wave operators W^{H^ + H^) on the waveguide with their 
analogues W^{Hy + Vr, i^r) on F: this is sketched in the last section. 

The method of approximating (zero-energy) resonant sequences 

We have a sequence of star-shaped waveguides {F^} ("^-fat" star graphs) which shrink 
to a star-shaped graph F as e — )■ 0. The method we describe allows the determination 
of the boundary conditions at the vertex of F, which are needed in order to define 
the effective dynamics on the graph that, for small £, approximates the low-energy 
dynamics generated by H^. These boundary conditions depend on the shape of F^ in 
the neighborhood of the vertex region Q^. The dependence is through the spectra of the 
auxiliary operators Hl^^ acting on the mesoscopic region F^^^^. 

We will see in the next paragraph that the control of the effective wavefunction ^ 
on the limit graph F near the vertex is possible if we consider a particular class of initial 
states (low-energy states) and choose wisely the parameters e, L. But once we restrict 
to this class of initial states, the gluing conditions at the vertex of the graph depend 
on the geometry of the waveguide and nothing else: the dependence is through the 
asymptotics (in the parameter L) of the eigenvalues of the rescaled mesoscopic region. 

The wavefunctions near the vertex 

Let ip'^ G L^(F^), smooth, normalized, energy-bounded state orthogonal to the bound 
states of F^. Consider the finite cylinders 

Bf^E'jX[0,i] j = l...n 

contained in the mesoscopic region: {£} x is part of the boundary of the mesoscopic 
region. Notice that (up to a zero-measure boundary) 

{UjBf) n = . 

The function \ B^'^ can be Fourier decomposed as usual: 
We can prove the following 

Theorem 2 Assume that the spectral gap condition holds. Ifip^ belongs to P(o,£;]iv^(F^), 
then the L°°-norm of 

= \[o/] 

is bounded by Ci ■ ^1'^ -|- C2 • £^/'^L~^'^~'^^f'^ , where Ci, C2 are numerical constants. 
♦ 

An important consequence of this theorem is the following. Let be the weak limit 
of any convergent subsequence {^1"}. It is in H^^'^(F\{V^}), therefore each component 
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j' = 1 . . . n is right-continuous at the vertex V . Since its approximants are uniformly 
bounded near the vertex by an infinitesimal quantity, it follows that 

l...n. 



lim = j 

a;-)-0+ 



Proof 



ijj^ is a state with energy bounded by E. Application of the spectral condition, 
together with (quasi)-orthogonality to the negative-energy states of H^^^^, tells us 
that 



e , /fo (e_\r 

£2 \l) 



^1 + 



yUO > 0, 7 > 



This inequality can be written in a more useful form: 

^ ^'^^ JTir,.t J Fir,., 



Now, we recall that 



f (ivv'"p-/x!iV''n >o 

^ Tout 



(this is easily seen by expanding ip^ in Fourier series) and so 

/ - f^wrf) < [ (ivv'i' - film') < E 

from which we deduce 

m' 



E 



This states that the norm of i/j Ivint (^'Hd hence of \Bj) is bounded by a suitable 
power of £ (not uniformly in E) . 

For any smooth function / on the interval [a, b] and any point xq G [o, b], the following 
Poincare-like estimate holds: 

f.6 



(b-a) \f{xo)\-\\Vfh.Vb^ 



< I l/C 

J a 



Consider the restriction of to ^: we write 



The remainder Rf ■ is the sum of the components of along the higher transverse 
modes. From the Poincare estimate applied to 

i\\^'{xq)\-\\V^'\\l2-1^/'^X < /V'l^< I W?-^ 



Since 



^ we recall that these states approximate the bound states of up to an exponentially small error, 
see the Definition of Hf^^. We omit the exponentially small errors from the formulas, to make them 
easier to read. 
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1 — I 

2 



2-7 • 
1\~ 



is finite, we obtain an estimate for sup l^'^l = ||*^||oo in terms of ^, e: 

W{.xo)\ < E ■ £1/2 + ^/E/JTo ■ 6^/2 (^^' 
Using e^/^ = ^^/^ • L~^/^ this last equation can be rewritten as: 

If 7 < 2, the term in curly brackets is bounded by E + ^E/ ii-y as £ ^ 0, irrespective 
of the particular choice of the parameter L that we have introduced. If 7 > 2 then 
we need to restrict the choice of L, in order to bound |^''(xo)| with an infinitesimal 
quantity. In fact, by choosing L = for a sufficiently close to 0, the last equation 
becomes 

|^^(xo)| < E ■ £i-" + ^/e/JTo ■ 

and the right hand side clearly vanishes in the limit e — t- for any value of 7, if a is 
suitably chosen. 

In the following it will be useful to know something about the restriction of ip'^ to 

the boundary of Vlf^{. 

Theorem 3 Assume that the spectral gap condition holds. Ifip^ belongs to P(o,B]L^(r'^), 
then the LF'-norm of 



is bounded by Ci ■ H}!"^ + C2 ■ £i/2_L-(2-7)/2^ where C\, C2 are numerical constants. The 
constants are not necessarily the same of theorem 2. 
♦ 

Proof 

Consider the sum of the higher modes ^'1, ^'f...: using in the Poincare-like 
estimate and summing up from m = 2 to 00, we get 

2 00 



m— 9 V J J m— 9 



||2 
■ m I 

m=2 ^ •> / m=2 

(all norms are L^-norms). We use the bounds obtained in Theorem 1 for the sum of 
the higher modes: the proof then mimics the proof of theorem 2. 

Another simple modification of Theorem 2 allows to state that, if the resonant state 
condition holds: 

Theorem 4 Assume that the resonant sequence condition holds. If ip^ belongs to 
P(o,E]L^(r''), define 



C.T/ — 
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and 

Then, the L°° norm of'^^, restricted to the segment [0,i], is bounded by Ci ■ E^^"^ + C2 ■ 
£i/22^-(2-7)/2^ yjliere Ci, C2 are numerical constants. Recall that {i^^j^+i) ^■^ resonant 
sequence ofT]:^^. 
♦ 

Pictorially, we may say that -0^ "relaxes" to V'fe+i- 
Dirichlet dynamics on graphs 

In the spectral gap case, our argument is already sufficient to estabhsh that the hmit 
dynamics on the star graph is the decoupled one with Dirichlet boundary conditions 
at the vertex. Let 4>'^{t) be a solution of the Schrodinger equation in P(o,E]L^(r^). We 
know that 

r{t) \^s,^^>\{t)®x\+Rm. 

out 

To our purposes, R\ can be neglected. The function ^\ has a weak hmit ^, for £ — > 0; the 
limit exists by compactness (choose a subsequence), and it satisfies the weak Schrodinger 
equation on r\{V^}. By Theorem 2, ^ satisfies Dirichlet boundary conditions at the 
vertex V: 

*^.(0) = j = l...n 

and this fixes the unitary evolution of ^ on the limit graph F completely. 

Since any convergent subsequence has the same limit at time zero, and the limit 
evolves with the same equation, then the whole family {^f } converges to the same limit 
^, at all times. 

Kirchhoff dynamics on graphs 

We have seen that the eigenstates of H^^^ which satisfy Equation 1 (the spectral gap 
condition) give in the limit e ^ a contribution to the limit wavefunction which is 
negligible near the vertex. The limit behaviour at the vertex other than Dirichlet 
conditions is entirely due to the resonant sequence of eigenstates; we must therefore 
analyze in detail this resonant sequence. In particular, we are interested in finding how 
the resonant sequence determines the parameters which characterize the self-adjoint 
extension of the Laplacian on the graph. 

The resonant wavefunction V'^+i Let the complex numbers (31'^ . . . (3^'^ be defined by 
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where the term rj stands for the sum of the higher modes x^, m > 2. We will choose 
the normalization of V'fc+i so that 

n 

With this normalization the numbers /3J'^ depend on e and i only through the ratio 
e/£ = 1/L: this is a consequence of the homotheticity of the family of waveguides {F-}. 
Define 

6^ = lim B^/ 7 = 1 . . . n : 
(we assume that the limit exists). Note that by continuity 



Consider the resonant eigenvalue A^'^^. Set a = e/i; a = 1/L, so that in our setting 
(T <^ 1. By a simple scaling argument, 

where Afe+i(o") is the {k + l)-th eigenvalue of the rescaled mesoscopic region of thickness 
e = 1 and length L = 1/a. Define 

^ Afc+i(a) - Afc+i(0) 

o--)-o a 

(we assume that this limit exists, too). 9 is (minus) the derivative of the resonant 
eigenvalue A^+i with respect to the parameter a. We point out that 9 is nonnegative. 

The following estimates are useful to recover the boundary conditions at the vertex. 
Consider the function ijj'^ defined on F^^^, solution of the inhomogeneous Dirichlet 
problem 

—Au — — 

t{£}xs|= • X! j = l...n 

Theorem 5 The following formula holds true: 

^ = hm £ • 9''^ 9'^^ = - I )*d^^^ . 



Proof 

Observe that V'fc+i is the solution to an inhomogeneous problem similar to the one 
solved by -0"^, but with different boundary data on {£} x S^: 

/i^^fcfi \{i}x^= (if ■X\ + ■■■ 

Since both the boundary data and the eigenvalues of the two systems converge, we 
expect the solutions to converge to each other in the limit £ — > in a suitable Sobolev 
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topology. If the boundary data converge in L^, converges to uniformly, but 

their derivatives (in particular, the normal derivatives at the boundary) in general 
do not. 

Consider the Gauss-Green identity 

Jn Ja JdQ. Jdci 

and apply it for / = (V'"^)*, g = i'l'+i^ ^ = ^int', using the equations 

we get 

[ Af.g- [ f.Ag = X^^, I 
t/ r2 •/ r2 J Tint 

which is equal to 

f dyf -g - i f -dyg 

JdQ JdCi 



Note that 



and 



J- int. *J 1 i.Tt.t 



lim:^^^ = V|/3,p = l. 



Thus 



lime-r-^ = lim-e-^- ( Xli.) 



k+1 ^k+l ) 



(note that A^'^ = 0, this term was introduced to emphasize the fact that the 
definition of 9 involves "derivative" of the resonant energy level A^.'^^^ . 



Kirchhoff graphs The argument runs exactly in the same way as in the Dirichlet case, 
except for establishing the vertex conditions, which is what we are going to do. Let us 
consider a generic wavefunction in P(^q^e]L'^{^^)', recall that 

out 

by evaluating this expression on {£} x we see that the right hand side equals 
(dots stand for the negligible higher modes). By Theorem 4 we see that 

\{i}xT,^. = C% ■ Ipl'l^ t{€}xE| + • ■ ■ 
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and since 

the comparison of the coefficients of the fundamental mode Xi gives 

n.(^) = 4-/3f + ••• 

which, for the hmit wavef unction ^, reads 

Now consider the restriction of ^'^ to the mesoscopic region F^^^^. We make use of 
the Gauss-Green identity: 

[ Af-g- [ f-Ag= (f dj-g-(f f ■ d.g , 
Jn Jn Jdn Jon 

substituting Tint for (■0'^)* for / and ip^ for g. 

The volume integrals vanishes in the hmit £ — >■ 0. In fact, using that —Aip"^ — jj^ip'^ 

we obtain 

f Af.g- I f-Ag^ ! {Ar)* ■ ^ - f {i^')* ■ A^ 

= I {^y[-A-^^\]r. 

It is clear that [—A — /xf JV'^ belongs to P(o,E]-L^(r^) as '4>^ does: by dividing the integral 
into the region VL^ and the cylinders, it is easy to see that this whole volume term is 
vanishing like 0{€), or faster. 

Now consider the surface integrals. By Theorems 3 and 4, we know that at the 
boundary of dVint ip'^ is proportional to ipl'^i, plus negligible (in L^-norm) terms: 

r ^ 4 • i^ii, ■ 

Moreover 

' dj-g-i f-d.g=I {d.^y-r-{ {i^'Y-d.r 

an Jan Jdri„t Jdri„t 



We can substitute the Fourier expansion for ^0^, and as usual the only nontrivial 
contribution comes from the fundamental mode : apart from negligible contributions, 
the surface integral equals 



1 " 
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Remark The e^^ factor above is rather surprising: notice that we can factor 
it out if we use the microscopic coordinate X {x = eX), so that another comes 
out of the derivative term. We want to stress that, despite that the use of X rather 
than X may look curious, our aim is, given [e is small, but finite), to determine 
uniquely n+ 1 numbers fii . . . [in, {oi 6 / e\) which are sufficient to completely describe 
the effective dynamics on the hmit graph T. The fact that these numbers may be 
normalized differently, or be covariant with respect to reparametrizations of ^'f, is 
physically irrelevant. 

If we choose to work with the macroscopic coordinate X, we can easily take the 
limit of the last equation for £ ^ 0, and we obtain 

n 

To get quantitative estimates of the error terms, we have to assume some concrete 
estimate on the convergence of the parameters ^J'^, 9^'^ to their limits: for example, 

- p. = OiL^"), J = 1 . . . n; e- 9''^ -9 = 0(1-") 

for some exponent k gives an error 0{e~^L~'^). Summarizing, we established that the 
limit wavefunction ^ satisfies 

e;=i/3;-^;(o)+^^-c,. = o 

which is sufficient to define the domain of the self-adjoint Hamiltonian Hr- 
5. (Weak) convergence of resolvents 

In this section we compare the resolvents of the waveguide Hamiltonian — — A — //f 
(at finite e) and the graph Hamiltonian H^. We will discuss expficitly the Kirchhoff case, 
so we will consider the situation of waveguides satisfying the resonant state condition; 
the same argument, with straightforward changes, applies to waveguides with spectral 
gap. 

We consider a lift (constructed below) of the self-adjoint extension Hr of the free 
Laplacian on r\{V^}. We will prove that the resolvent {H^ — weakly converges to 
(iff ~ ^)~^ if s-iid oiily if the parameters which define Hr coincide with those determined 
by the resonant sequence. Convergence holds true for finite-energy states orthogonal to 
the bound states of F^. 

Notice that the set of self-adjoint extensions of the free Laplacian on the graph is 
much larger than the set of extensions that we obtain: our construction indicates that 
only this small sub-class can be obtained as a limit of the Dirichlet Laplacian on F^. 
We shall come back to this point in the conclusions. 
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According to our construction, the Hamiltonian is zero on states which have 
no component in the first transversal mode of the channels. Since these states belong 
to the continuous spectrum of H^, no stronger convergence is to be expected for the 
resolvents. 

Lifting graph Hamiltonians to 

Let us consider the self-adjoint Hamiltonian Hy on the graph F, extending the free 
Laplacian on the edges; the domain of the Hamiltonian is characterized by the set of 
numbers Pi . . . (3ni 9 va. the following way 



Given a function \E' : F — )■ C, in the domain 'D[H), we can lift it to a function on F^ 
with the operator J as follows: 



In this way, the whole domain V can be hfted to F^: on the closure of this Hilbert 
subspace, we may define the lifted Hamiltonian Hy- 



and on the complement of JV we simply set equal to zero. This means that, for 
instance, iff is null on all "high-energy" states ^ ® xtm > 2. 

A useful computation 

The following computation will be needed in the comparison of resolvents; for easier 
reference, we prove it separately and state the result at the end of the paragraph. Let (p 
belong to the domain of and t/j belong to the domain of iff, that is to say, tp — J^, 
for ^ e V{Hy). We want to compute the difference 





on f; 



Let us compute with {H^(p, ip). Since vanishes on F^^^^ this reduces to an integral 



out' 




A simple integration by parts leads to 





We write (/? = $ (8) xf + . . . (dots stand for higher transverse modes) : then 
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Now, we compute ((/?, i^f ■0) ■ Again, this is an integral on F^^j only: 

Substituting the Fourier decomposition of </? and integrating by parts, we get 

n „ 

(v.,i/ft^) = ^$*(o).vi/;.(o)+ [^'ym' . 

7=1 ''^ 



Recall that $j(0) = c$ ■ = c# ■ 13 j + 0{e the sum above can be simphfied as 

follows 



i=i \j=i / 

It is immediate to verify that equals Le,e (V</7*V'0 - A*i</7*'0); so we 

conclude that 

(i/V,^) - {if.HI.^) = . 1^/3* . $;.(0) + ^ ■ + 0(£-iL-'^) ; 

the term in curly brackets vanishes (like 0{l^/'^)-\-0{e~^L~'^)) if and only if the boundary 
conditions of Hy at the vertex are chosen as those implied by the resonant sequence. 

The resolvents of and 

We are now in position to compute the difference of the resolvents of and These 
operators are bounded, so that it suffices to control their difference on a suitable dense 
subset of the Hilbert space L^(F^). 

Let ^ be a complex number with Jm(z) 7^ 0; let be a fixed positive number. Let 
us consider (/?,-?/'€ M^^'^(r^), normalized, with ip in the domain of /Jp. Both ip and ip 
have energy less than and are orthogonal to the bound states of F*^ (i.e. they have 
components in the continuous spectrum only). We prove the following 

Theorem 6 Consider the expectation value 

F'(^.V.)-(^.(^-^)* 

where (p, ip are normalized states in L^(F'^), with ip orthogonal to the bound states and 
of finite energy: then 

♦ 



We remark that the constants multiplying the various powers of e and I actually 
depend on the energy and blow up when E ^ 00. 

Proof 
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We will start with the formula for the difference of the resolvents: 

1 -m-H-y ' 



-z H^-z -z^ ' 'Hf.-z 

by substituting this formula into our expectation value, we get 

where 

= - z)-\ V. = m - z)-^^ . 

We remind that in the computation of the resolvent formula, iJp is meant to be 
applied to V'z, while H'^ is applied to ipz'- 



We will compute this formula for if G T>{H^) and for = J^, * G T){Hr). We will 
briefly sketch how to treat the case of ip supported in T-^^ and supported in F^^j 
with high-energy components: the theorem then follows by density. 
Notice that = J^z, where *2 = (Hr - z)~^^. Clearly, G T^{Ht). Moreover, 
ipz G V{H^). Thus, using the formula that we deduced in the previous paragraph. 

If ^ is compactly supported in T^^^, the action of on ^ is trivial. The remaining 
integral on F^-^^ is estimated as usual, and is negligible. If ip is an high-energy state, 
the action of Hf, is again trivial, and we notice that the components of tp along the 
higher transverse modes are suppressed: hence, the integral vanishes as £ — >■ and 

This theorem implies that on scattering states with finite energy the difference of 
the resolvents goes to zero if £ — > and £ — > 0. 



The scattering problem 

Consider a waveguide in R^, and a smooth bounded fast decaying potential V defined 
on IR^. Let as usual denote by the operator — A — //f , where A is the Laplacian with 
Dirichlet boundary conditions on dV^. Assume that admits exactly one resonant 
sequence and let — — Ap, where Ar is the Laplacian on F with vertex conditions 
adapted to the resonant sequence. Both H'^ and Hr have the positive real axis as 
absolutely continuous spectrum. In addition, may have bound states. One can 
define the wave operators 

where Pa.c. is the projection onto the absolutely continuous spectrum, and 
W^(Hr + Vr, Hr) = s - lim e^t(Hr+Vr)^-itHr 

where Vr is the evaluation of V on the graph F. Asymptotic completeness holds in 
both cases, and one can define unitary S- matrices S'^{V), Sr{V). It is easy to find the 
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explicit form of W^{Hy + Vr,-ffr) with perturbative methods of stationary scattering 
theory, since the generahzed eigenfunctions of are known. It is more difficult to give 
an explicit form of W^{H^ + V,H^)] therefore it is useful to have an approximation 
of it in terms of W^{Hy + Vr, -f^r)) with error estimates. We sketch here the way this 
approximation is obtained by the method of resonant sequences: it is essentially an 
adaptation of the chain rule. Details will be given elsewhere. 

We start with the identity, vahd for all t and all positive values of the small 
parameter e: 

^-itHI,^t{HI.+V)^-it{HI.+V)^it{H'+V) ^ ^-itHl.^tH'^-itH'^it(H'+V) 

(clearly, we use the lift of Hr; the wave operators for Hf^ arc the obvious extension 
of the wave operators of Hr). By Wiener's lemma, for every E > one has 

hm hm -e-^*^f e^*^^n£ = Ue hm hm ie-^*(^f +^)e^*(^^+^)nE = 
£^0 T-^oo T £->o r->oo T 

where He is the projection operator on states in which < E. Since strong limits can 

be composed one has for any E 

WmllEW^iH^^ + V,H^^)IIe = UE\imW^(m + V, m)UE . 
It is easy to recognize in the last term W^lHr+Vr, Hr), after the natural identifications. 

6. Comparison with previous results 

We compare briefly our approach and results with previous work on the same problem. 
There are many published papers on the subject, here we limit ourselves to only some 
of them. 

As we remarked, our main emphasis is on the fact that the limit dynamics at a 
vertex of the graph depends on the spectral properties of a sequence of Schrodinger 
operators defined on neighborhoods of the vertex in the fattened graph, and that in 
order to obtain in the limit boundary conditions other than Dirichlet it is necessary 
that a resonating sequence exists for this sequence of operators. 

The role of some sort of resonance has been stressed by B. Pavlov (see (Harmer 
et al. 2007) and previous papers). In (Albeverio et al. 2007) this is rigorously proved 
in the case of a graph shrinking to a line with a sharp bending, in the context of limit 
point interactions. In this context the asymptotes at ±oo of a zero-energy resonance 
provide, in the suitable scaling limit, a connection between the two sides of the singular 
point on the line. This role of the resonating sequence also in the case of graphs has 
been our guiding idea. 

In a very interesting paper (Post 2005) O. Post has proved (among other things) 
that some shrinking of one of the tubes towards the vertex to which it is attached is 
enough to provide Dirichlet boundary conditions in the limit for the edge corresponding 
to that tube. Our results imply that the shrinking is not necessary if a resonant sequence 
does not exist. 
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In another relevant paper S. Molchanov and B. Vainberg (Molchanov & Vainberg 
2007) prove (among other things) continuity of the scattering matrix under the shrinking 
operation e — )■ for "fattened graphs" . This is also sufficient to establish a weak form of 
resolvent convergence, very similar to the one that we propose. This abstract result does 
not give information on the role of the geometry of a neighborhood of the vertex and 
therefore does not use information that comes from the specific geometries suggested by 
physical examples, e.g. from the density of conduction electrons on aromatic molecules. 

A refined analysis of the structure of eigenfunctions can be found in the very 
interesting and relevant paper of D. Grieser (Grieser 2008a), which again exploits the 
strategy of studying the scattering matrix. The paper of Grieser is very detailed and 
full of relevant results. In particular one can extract from it that in the case of "fattened 
graphs" which shrink uniformly in e a resonant sequence has energy exponentially close 
to the continuum treshold. While our approach is local and relies on energy esitmates 
and Neumann-Robin bracketing, the approach of (Grieser 2008a) relies on connecting 
smoothly the eigenfunction of the internal region with scattering waves. In spite of the 
different estimates, there are points in common between the two approaches. It would in 
particular be interesting to extend Grieser's results to sets with nonuniform scaling, 
which have a geometry more in line with physical data. 

Note that a suitable notion of resolvent convergence can be achieved both within 
our approach, and through the continuity in e of scattering matrices. We think that 
our approach makes the role of the resonant state very transparent: we introduce 
the resonant sequence, which is a kind of "local" (restricted to a neighbourhood of 
the junction) notion of resonance, and we identify the parameters of the limit graph 
Hamiltonian in terms of the asymptotic values and the energy of the resonant sequence. 
Moreover, the dependence of the limit Hamiltonian on the geometry of the junction 
is actually through the spectrum of the mesoscopic region, which is a significant 
simplification. 

In our paper we give an estimate of the error one makes as a function of e on the 
evaluation of relevant physical quantities when using the limit dynamics; this can be 
done in all approaches which take into account the geometric shape of the neighborhoods. 

7. Conclusions and generalizations 

Let be an £-thin, star-shaped waveguide and consider a Schrodinger evolution 
generated by the Hamiltonian — —A — //^, where A is the Dirichlet Laplacian 
on r^; choose /xf so that the onset of the continuous spectrum is 0. Define a sequence of 
auxiliary (disconnected) waveguides F^^^, isometric to F x E^, together with a sequence of 
lifted Hamiltonians under the assumption that there exist (at most one) resonant 
sequence of eigenstates, we prove that the resolvent of Hf. converges weakly to the 
resolvent of H^, on low-energy states. Estimates of the difference are given, as well 
as estimates for the approximation of observables such as density and (longitudinal) 
current. 
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The argument is constructive: the parameters that identify are uniquely 
determined in terms of simple properties of the resonant sequence, which in turn depends 
on the geometry of the vertex. When the same results hold for free Hamiltonians 
perturbed by an external potential V ^ the wave operators W^{H^ + V, H^), "projected" 
on low-energy states, indeed reduce to their analogues W^{Hy + Vr, Hy) on the graph 
r, and this provides a useful approximation formula. If the resonant sequence is absent, 
we prove that the limit Quantum Graph is the decoupled one. The case with more than 
one resonant sequence is not addressed here. 

We point out that the effective dynamics on metric graphs that we obtain with our 
approximation method belongs to a subset of all possible self-adjoint extensions of the 
free Laplacian on r\{\^} that are characterized by the boundary conditions 

n 

where j3j is the asymptotic value of the resonance in the j-th branch of the waveguide 
and 9 is the derivative of the resonant energy. In view of applying our methods to study 
diffusion on thin waveguides (Preidlin & Wentzell 1993), we remark that these, whenever 
the can be chosen nonncgative, are precisely the conditions under which generates 
a positivity preserving contraction semigroup on F (with continuous trajectories when 
Pj = Pk) (Kostrykin et al. 2008). If P^o = for some the corresponding branch carries 
an independent process. On the other hand if H"^ has no bound states it generates 
a positivity preserving contraction semigroup on F^; it is natural to expect that the 
stochastic process generated by when averaged in the transversal direction converges 
to the process generated by Hr when £ — > 0. We shall come back to this issue in a future 
publication. 

We conclude with a few remarks. The result described here holds also in the case 
in which the arms of the waveguide are e-neighbourhoods of smooth curves, provided 
that their sections are compact, of linear size in e and the bottom of their Dirichlet 
spectrum is nondegenerate. The method can also be used, with simple modifications, 
to treat generic waveguides which are ^-neighbourhoods of graphs with a locally finite 
number of vertices (better: the distance between any couple of vertices is bounded from 
below by a positive constant) (Costa 2010). 

Even if the method described here is used for a sequence of homotetic thin 
waveguides, it can be applied to the case where the shape of the £-waveguide depends 
on the parameter e. In particular, it can be applied to study a two-dimensional strip 
with a bend which becomes sharp as e — )■ 0. In (Albeverio et al. 2007) this problem 
was studied in two dimensions, the curved strip was shrunk to a smooth curve and then 
the curve was sharply bent. In this case the motion on the curve is described by the 
free Schrodinger equation plus an attractive potential (due to the curvature) that in the 
limit £ — > gets deeper and concentrates around the origin (the sharp bending point). 
If the parameters are conveniently chosen one can obtain a zero-energy resonance and 
a limit motion characterized by Kirchhoff-like boundary conditions, with parameters 
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depending on the resonance. For all other choices one gets decouphng conditions at the 
origin. 

As for concrete examples of resonant sequence, taking inspiration from pictures of 
aromatic molecules, we will show that one can have a resonant sequence for a waveguide 
composed of a thick spherical shell Re < \x\ < R'^ attached smoothly to cylinders of 
base where the parameters are chosen in such a way that the first eigenvalue of the 
Dirichlet Laplacian in the spherical shell is equal to the first eigenvalue of the Dirichlet 
Laplacian on E^. 
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